Phonon-mediated anisotropic superconductivity in the Y and Lu nickel borocarbides 
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We present scanning tunneling spectroscopy and microscopy measurements at low temperatures 
in the borocarbide materials RNi2B2C (R=Y, Lu). The characteristic strong coupling structure due 
to the pairing interaction is unambiguously resolved in the superconducting density of states. It is 
located at the superconducting gap plus the energy corresponding to a phonon mode identified in 
previous neutron scattering experiments. These measurements also show that this mode is coupled 
to the electrons through a highly anisotropic electron-phonon interaction originated by a nesting 
feature of the Fermi surface. Our experiments, from which we can extract a large electron-phonon 
coupling parameter A (between 0.5 and 0.8), demonstrate that this anisotropic electron-phonon 
coupling has an essential contribution to the pairing interaction. The tunneling spectra show an 
anisotropic s-wave superconducting gap function. 

PACS numbers; 74.70.Dd, 74.62. Dh, 74.20.Mn 



In most known superconductors the formation of 
Cooper pairs is based on an attractive interaction me- 
diated by phonons. However, Cooper pairing driven by 
other bosonic excitations has also attracted much atten- 
tion due to its fundamental interest. Some of the pro- 
posed mechanisms have been used to interpret the su- 
perconducting behavior of different compounds discov- 
ered during the last decades. But clear microscopic in- 
formation is very much needed in this field. The situ- 
ation is especially puzzUng in the borocarbide materi- 
als (RNiaBaC, R=Y, Lu, Tm, Er, Ho, Dy), where re- 
cent experiments have shown that the question about 
the pairing interaction, initially thought to be conven- 
tional electron-phonon coupling||l|, ^, ^, is far from be- 
ing understood. These compounds show moderate criti- 
cal temperatures (between 6K and 16.5K[^ |5|) and very 
interesting phase diagrams where superconductivity co- 
exists with antiferromagnetic order (when R is a mag- 
netic rare earth, R=Tm, Er, Ho, Dy). The behaviors ob- 
served in the thermal conductivity |^, photoemission 
spectroscopy!^ specific heat||, |l^, microwave surface 
impedance P, and Raman-scattering [p2| experiments 
of the non-magnetic Y and Lu borocarbides (which also 
present the highest critical temperatures of 15.5 K and 
16.5 K respectively) show that the superconducting gap 
is highly anisotropic. Indirectly, this could be related to 
an also anisotropic pairing interaction, but no experiment 
has given an indication of its nature. 

On the other hand, Fermi surface nesting seems to be a 
general feature of the whole family of borocarbide mate- 
rials. It produces Kohn anomalies and has directly been 
observed using angular correlation of electron-positron 
annihilation radiation in YNi2B2C and LuNi2B2C|13, 



p^ , [T^ . Moreover, the antiferromagnetic ordering found 
in the borocarbides with a magnetic rare earth has a 
wavevector very close to the nesting vector Q« (0.5, 0, 0) 
|l6| , [T^ (in R=Er, Ho and Dy; in Tm it appears in a 
magnetic field phase). Nesting is generally considered to 
play against the formation of a phonon mediated super- 
conducting state, as evidenced in other families of super- 
conductors (e.g. in the dichalchogenides ||l8|, |l9[). Here 
we present measurements of tunneling spectroscopy of 
the non-magnetic YNi2B2C and LuNi2B2C, which show 
that the highly anisotropic electron-phonon interaction 
produced by Fermi surface nesting drives the system to 
an also anisotropic superconducting state. We give new 
insight into the nature of the pairing mechanism, the or- 
der parameter symmetry and the gap anisotropy. 

Tunneling spectroscopy is in principle one of the most 
powerful experimental tools to investigate the anisotropy 
of the superconducting gap and to obtain information 
about the pairing mechanism|Q. We use the same scan- 
ning tunneling microscopy and spectroscopy (STM/STS) 
set-up as in ref . ||2l| , where we studied the magnetic 
TmNi2B2C, with an improved resolution allowing now 
measurements down to 0.5K. We have also characterized 
other materials with the same set-up (Al[^, Pb and 
NbSe2) and demonstrated that this is indeed the actual 
temperature of tip and sample and that we do not need 
any additional pair breaking parameter to explain 
our data. Previous works about the spectroscopy of boro- 
carbides had a lower resolution in energy |2^, p5[ . 
The sample is broken in air on the sample holder of the 
STM and cooled down immediately. The resulting sur- 
face presents the same topology as for the Tm borocar- 
bide, which consists of inclined planes and bumps, typi- 
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cal of a conchoidal fracture, presenting no clear crystallo- 
graphic orientation[^. So it is crucial to characterize the 
superconducting behavior in well-differentiated regions of 
the sample. We use a home made x-y table that gives 
the possibility to change in-situ the scanning window (of 
Ixl/im^) in an area of 2x2mm^. We measured three 
different samples (in three different cool-downs) of each 
compound. In the case of LuNi2B2C, all of them were 
freshly broken pieces of the same single crystal grown by 
a flux technique described in ref. Q. The YNi2B2C sam- 
ples came from two different single crystals, one grown 
by the same flux technique and the other in an image fur- 
nace. The tunneling conditions were always very good, 
with high valued measured work functions of several eV. 
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FIG. 1: a) Tunneling conductance spectra (taken at 0.5K 
and normalized to the conductance value at high bias volt- 
age) where electron-phonon coupling features at high voltages 
clearly appear, b) Histogram counting the number of times 
we find such a feature (within a sample of a hundred spectra), 
with the maximum value of the derivative of the conductance 
at the voltage reported on the x axis. 

The tunneling differential conductance, a—dl/dV, be- 
tween a normal metal and a superconductor gives a di- 
rect (temperature smeared) measurement of the super- 
conducting density of states. Therefore, data show the 
eventual presence of very low energy excitations within 
the superconducting gap, 2A, and a high quasiparticle 
peak at voltages close to the gap. But if the measure- 
ment is sufficiently precise, one can also try to resolve 
tiny features at voltages corresponding to the sum of A 
and the characteristic energy of the bosonic excitations 
leading to superconductivity, in order to obtain informa- 
tion about the nature of these excitations] 20, 2l, p7| 



We 

could indeed observe these features in the tunneling con- 
ductance spectra for the Y and Lu borocarbides. In Fig.l 
we show a typical curve together with a histogram where 



the number of times that such features are observed is 
presented as a function of the bias voltage (their voltage 
position is determined from the maximum of the deriva- 
tive of the conductance above the gap[^ ^ |27|). Most 
of them appear centered at a voltage between 6 and 7mV. 
Substracting the values of A, as estimated below, the 
maxima become centered between 4 and 5 mV. This is 
precisely the energy at which a high peak develops in neu- 
tron scattering experiments due to a low energy phonon 
mode|l^, An estimation of the electron phonon cou- 
pling constant A from our data gives rather high val- 
ues between 0.5 and 0.8 in both compounds [p8[. On the 
other hand, previous estimations using thermodynamic 
measurements, which take into account the Eliashberg 
function at all relevant phonon frequencies, but which do 
not directly specify the nature or spectral weight of each 
mode, have given values of A between 0.75 and 1.2(29[ |o). 
The fact that we find values of A smaller but very close to 
the ones estimated with thermodynamic measurements 
means that the structures in the Eliashberg function as- 
sociated with higher energy phonon modes Q are 
more spread-out in energy than the peak corresponding 
to the 4mV phonon mode and fall below our experimen- 
tal resolution. It also means that the low energy phonon 
mode measured here is essential in the formation of su- 
perconducting correlations. 

The most striking point however is that the phononic 
density of states at these energies results from a mode 
having a wavevector comparable to the nesting vector Q 
(0.5,0,0) |l7). Neutron scattering experiments show 
indeed pronounced Kohn anomalies [Esl p^ , where the 
nesting feature of the Fermi surfacefla, |17[ leads to a 
significant softening, when decreasing temperature, of 
the low-lying transverse phonon branches at wavevectors 
close to (0.5,0,0). This behavior results in a strong and 
highly anisotropic electron-phonon coupling. Our experi- 
ment demonstrates for the first time that this anisotropic 
electron-phonon interaction, produced by Fermi surface 
nesting, leads to superconducting correlations, having an 
important contribution to the total electron-phonon cou- 
pling constant A. Correspondingly, we can expect the 
superconducting gap to be also highly anisotropic. 

Tunneling conductance measurements done with a 
STM give the superconducting density of states averaged 
over part (but not all) of the Fermi surface, depending 
on the relative position of the tip onto the sample ||3^. In 
most cases we measure tunneling conductance curves as 
the ones shown in Fig. 2, where the conductance is zero 
below O.SmV and then increases up to a peak located at 
2.6mV in Lu and 2.3mV in Y. These curves cannot be 
fitted by conventional BCS theory. The discrepancy with 
most simple BCS theory is not due to a lifetime smeared 
BCS density of states, which leads to a non zero den- 
sity of states at low energies, not observed in our data. 
The situation is clearly very different from the one found 
in the very similar borocarbide material TmNi2B2C [Ell, 
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FIG. 2: Tunneling conductance spectra in the Lu (a) and Y 
(b) borocarbides at 0.5K normalized to the conductance value 
at high bias voltage (tunneling resistance of 2MQ). The full 
lines are fits to conventional BCS theory (A=2.3meV). 



where the spectra can be fitted by BCS theory. The form 
of the curves in Fig. 2 shows that the electrons contribut- 
ing to the tunnehng current at this precise location come 
from parts of the Fermi surface with a continuous dis- 
tribution of values of the superconducting gap|^, |5). 
This reveals that the gap function must be anisotropic. 
Moreover, the spectra in (Y,Lu)Ni2B2C remain with the 
same shape over regions much larger than the coherence 
length (usually dimensions about 200 x 200nm^ or bigger 
whereas ^'--^Tnm Q), maintaining their overall form when 
we increase the temperature, and becoming completely 
flat at the bulk critical temperature. Similar observa- 
tions have been made in the anisotropic superconductor 
NbSe2 where the spectra have, qualitatively, the 

same form. 

On the other hand, although the behavior shown in 
Fig. 2 appears most frequently on the surface, we can also 
find different topographical regions of similar sizes with 
very different spectra, as shown in Fig. 3a. In Fig. 3b we 
compare the temperature dependence of A, as estimated 
from the voltage position of the maximum of dcr/dV, and 
in Fig. 3c the corresponding tunneling density of states at 
0.5K in three different locations of YNi2B2C (the same 
is observed in LuNi2B2C). Clearly, the differences in the 
spectra are associated with smaller values of the super- 
conducting gap and also a reduced critical temperature. 
Remarkably, the mean value of the superconducting gaps 
shown in Fig. 3b follows well the temperature dependence 
predicted by simple BCS theory (lines). What is more, 
the curves corresponding to a smaller critical temper- 
ature in Fig. 3c (top curve) approach much better an 
isotropic BCS s-wave behavior than those correspond- 
ing to the bulk Tc (bottom curve) . This must be related 
to changes in the anisotropy as a function of the local 



Tc. To obtain more quantitative information about this 
effect we have fitted the experiment to a modified BCS 
s-wave density of states assuming that the dispersion in 
the values of the superconducting gap can be modelled by 
a gaussian distribution centered around Aq with a width 
of e (a similar approach has been previously used in the 
anisotropic superconductor NbSe2[Q). The dashed lines 
in Fig. 3c give the best fit to this model. The agreement 
with the experiment is much better in regions with a 
smaller Tc, where e decreases. The inset in Fig. 3b repre- 
sents the dependence of the estimated anisotropy, e/Ag, 
as a function of the measured critical temperature in sev- 
eral topographical regions of the same samples. Note that 
these values do not give an indication of the whole distri- 
bution of values of the superconducting gap, as a given 
STM spectrum is a local measurement which probes only 
part of the Fermi surface|Q. Therefore, our data are in 
good agreement to previous macroscopic measurements 

10,1,1,0, ini, 111, 
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FIG. 3: a) Topografic (left) and STS (right) images at 0.5K 
(YNi2B2C, aiOxSlOnm^). The latter figure is indicative of 
the superconducting density of states as a function of the posi- 
tion of the tip on the surface [j^ ] . The contrast is adjusted 
such that bright spots correspond to (t(V) similar to Fig. 2, 
and dark spots to cr(V) similar to the top curve in Fig. 3c. b) 
Temperature dependence of the mean value of the supercon- 
ducting gap (points) in YNi2B2C. We take A as the voltage 
where da(V)/dV is maximum. Lines give the temperature 
dependence of the superconducting gap expected in BCS the- 
ory. The inset represent the variation of e/Ao (see text) with 
Tc for YNi2B2C (circles) and LuNi2B2C (crosses), c) Spectra 
taken at 0.5K for each temperature run (same color for the 
straight lines as points in b), shifted for clarity. Dashed lines 
are fits to the BCS theory with anisotropy (see text), where, 
from top to bottom, Ao=l.lmV and e=0.2mV, Ao=1.3mV 
and e=0.32mV and Ao=1.8mV and e=0.5mV. 

Anisotropic superconductivity is expected to be very 
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sensitive to even non-magnetic defects that reduce Tc 
and decrease the anisotropy|^, 0. Therefore, in an 
anisotropic superconductor, local measurements in topo- 
graphically different positions on the surface can in prin- 
ciple show different forms of the superconducting gap and 
values of the local critical temperature, associated with 
the presence of defects lying near the surface or with an 
irregular topography[^. Isotropic superconductors, by 
contrast, present BCS spectra and the bulk critical tem- 
perature over the whole surface |2l|, |3^ . 

In (Y,Lu)Ni2B2C we observe a gradual decrease of the 
local critical temperature down to about half its bulk 
value and of the anisotropy (inset of Fig. 3b) in differ- 
ent locations which also show a different topography 
(Fig. 3a). Our measurements show a direct correlation 
between the local depression of Tc with a local change 
of the form of the gap. This also demonstrates that the 
superconducting gap must be highly anisotropic, and it 
provides an additional test of the nature of the order 
parameter in these materials. If the anisotropy is due 
to a d-wave order parameter, the defects tend to sup- 
press superconductivity altogether and low energy ex- 
citations appear, filling the density of states below the 
gap. If, however, the superconducting wave function 
is anisotropic s-wave, the defects tend to suppress the 
anisotropy, leading to a more isotropic gap and a de- 
creased critical temperature Eq, 41 . The curves plotted 



in Fig. 3c, with no low energy excitations, definitely im- 
ply that superconductivity in YNi2B2C and LuNi2B2C 
is highly anisotropic but s-wave. 

Using the same experimental protocol, we find a com- 
pletely different behavior in the chemically very similar 
magnetic TmNi2B2C compound (Tc^llK), where BCS 
s-wave like spectra without a significant anisotropy (and 
bulk critical temperature) are measured over the whole 
surface]^. However, we can state that the surface is 
analogous in all three compounds, because the measured 
work function and the topographical images are similar. 
This could be interpreted as an indication that mecha- 
nisms inducing anisotropy in YNi2B2C and LuNi2B2C 
are not operating in this borocarbide, or that the intrin- 
sic magnetic disorder already smears out homogeneously 
the anisotropy. 

In conclusion, we have studied the tunneling spec- 
troscopy in the non-magnetic borocarbides using high 
resolution STM/STS. We have been able to character- 
ize important microscopic aspects of the superconduct- 
ing state, which is an anisotropic s-wave state, where a 
significant part of the electron-phonon coupling leading 
to superconducting correlations is also highly anisotropic 
and due to soft phonons. The demonstration of this new 
mechanism strengthens the hope for further discoveries 
in the area of new high-Tc superconducting materials. 
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